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ABSTRACT 

We present the largest and most detailed survey to date of the stellar populations in the outskirts of 
M31 based on the homogeneous analysis of 14 deep HST/ACS pointings spanning the range 11.5 kpc 
< Rproj P$ 45 kpc. Many of these pointings sample coherent substructure discovered in the course 
of the INT/WFC imaging survey of M31 while others sample more diffuse structure in the extended 
disk. We conduct a quantitative comparison of the resolved stellar populations in these fields and 
identify several striking trends. The color-magnitude diagrams (CMDs), which reach > 3 magnitudes 
below the red clump with high completeness, can be classified into two main categories based on 
their morphologies. 'Stream-like' fields, so named for their similarity to the CMD of the giant stellar 
stream, are characterized by a red clump that slants bluewards at fainter magnitudes and an extended 
horizontal branch. They show no evidence for young populations. On the other hand, 'disk-like' fields 
exhibit rounder red clumps with significant luminosity width, lack an obvious horizontal branch and 
show evidence for recent star formation (~ 0.25 — 2 Gyr ago). We compare the spatial and line-of-sight 
distribution of stream-like fields with a recent simulation of the giant stream progenitor orbit and find 
an excellent agreement. These fields are found across much of the inner halo of M31, and attest to the 
high degree of pollution caused by this event. Disk- like material resides in the extended disk structure 
of M31 and is detected here up to R pro j ~ 44 kpc; the uniform populations in these fields, including 
the ubiquitous presence of young populations, and the strong rotation reported elsewhere are most 
consistent with a scenario in which this structure has formed through heating and disruption of the 
existing thin disk, perhaps due to the impact of the giant stream progenitor. Our comparative analysis 
sheds new light on the likely composition of two of the ultra-deep pointings formerly presented as pure 
outer disk and pure halo in the literature. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: halo — galaxies: individual 
(M31) — galaxies: stellar content — galaxies: structure 



1. INTRODUCTION 

A key goal of modern astrophysics is to understand 
the formation history of galaxies like our own Milky 
Way. The favored paradigm of hierarchical assembly 
within a cold dark matter (ACDM) dominated uni- 
verse predicts that spheroidal galaxy components form 
through a repetitive process of galaxy mergers and 
the accretion of smaller subsystems w hile disks arise 
from t he smooth accretion of gas (e.g., IWhite fc Frenkl 
119911: ISteinmetz fc Navarrol I2002D L Recently, much 
ground has been gained in understanding the nature 
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of stellar halos built up by tidal stripping of accreted 
sub- halos through N-bqdy pl u s semi-ana l ytical mod- 
els ([Bullock fc Johnston! 120051: iFont et alj l2006aft a nd 
through numerical simulations (|Helmi fc Whitdll999[ ). 

Assuming the accreted systems possess a stellar com- 
ponent, the most visible evidence of this hierarchical for- 
mation is expected in the form of tidal tails and partially 
digested satellites in the e xtended halos of galaxies (e.g. 
iBullock fc John ston 2005). Theoretical work on the tidal 
destruction of satellites by their massive hosts has shown 
that the mixing time of the debris depends on the satel- 
lites mass and orbit. This timescale can be relatively 
short for debris in t he inner galaxy but m any gigayears 
for the outer galaxy ([Johnston et al.|[T996[ ). Searches for 
spatially coherent tidal features in the extended halos of 
galaxies thus provide a means to probe ACDM predic- 
tions on scales of individual galaxies. Once tidal features 
are found, detailed follow-up studies of their stellar con- 
tent and kinematics provide important constraints on the 
nature and number of objects which have merged. The 
Milky Way hosts at least one ma jor tidal stream, orig i- 
nating from the Sagittarius dwarf dMaiewski et al.12 003') , 
with many others suggested (e.g.. iGrillmair fc Dionatoa 
[20061 [Beiokurov et al.ll20"ol . 

It can be argued, however, that the best current lab- 
oratory for studying faint stellar substructure around 
galaxies is provided by our nearest giant neighbor, M31. 
Our external view of this system removes complications 
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Fig. 1. — Map showing the distribution of RG B stars around M31 
from the INT/WFC survey Hrwin et al.l l2005). The image spans 
~ 95 kpc X 125 kpc. The locations of our HST/ACS pointings are 
overlaid in red, each covering 0.8 kpc X 0.8 kpc at the distance of 
M31. 

due to line-of-sight and extinction/crowding effects that 
plague Milky Way studies and as a result the struc- 
tures observed in M31 are considerably easier to inter- 
pret. M31 has the been the subject of several deep wide- 
field ground-based surveys in the last decade and abun- 
dant low-surface brightness stellar substructure has been 
found. These surveys have the depth to resolve individ- 
ual red giant branch (RGB) stars in the halo of M31 and 
the analysis of their large-scale spatial density distribu- 
tion has allowed unprecedented surface brightness levels 
to be reached (Sy > 30 magnitudes per square arcsec- 
ond). The Isaac Newton Telescope Wide-Field Camera 
(INT/WFC) panoramic survey was the first of these sur- 
veys to be conducted and imaged the outer disk and in- 
ner halo of M31 out to a radius of ~ 55 kpc revealing 
copious substructure and tantalizing evidence fo r meta l- 
licity variation s (Figu re 1 , see also llbata et alJ (|2001h ; 
IFerguson et~aT1 (|2002ft and Hrwin etai] (|2005h ). In the 
southern quadrant, a giant stream (hereafter GS) was 
discovered falling into the center of M31 from behind 
(Hbata et al.ll2001b iMcConnachie et~aT1l2003l : llbata et all 
2004). Prominent stellar overdensities were also discov- 
ered along both major axes, beyond the extent of the 
bright disk. More recently, the entire southern quadrant 
of M31 has been imaged out to a radius of ~ 150 kpc, 
with an extension to the galaxy M33 at ~ 200 kpc, us- 
ing Mega cam on the Canada France Hawaii T elescope 
(CFHT) ([Martin et alJ [200(1 llbata etail [200l . This 
survey has revealed an additional complex system of even 
fainter tidal streams in the far outer halo, including a se- 



ries of azimuthal streams along the minor axis out to 
~ 120 kpc. The Megacam survey has probed the GS out 
to a projected radius of ~ 100 kpc and uncovered evi- 
dence for internal stellar populations variations. Current 
simulations suggest it can be explained as the trailing de- 
bris arm of a disrupted dwarf galaxy of mass ~ 10 9 M Q 
whose leading arm m ay have wrapped around the inner 
ralaxv at least twice ([Mori fc Richl (|2008fl ; iFardal et ail 
2007) and references therein). The highly radial orbit of 
the s tream suggests it passes very close to the center of 
M31 ([Ibata et al.ll2004f ) and hence almost certainly inter- 
acts with the disk. This may represent one of the most 
significant accretion events since the initial formation of 
M31, providing us a fortuitous glimpse of ongoing mass 
assembly in the z = universe. The presence of the GS 
has raised many questions: What role has this particular 
event played in the late formation history of M31? How 
much of the other inner halo debris has come from the 
stream? What damage has the stream progenitor done to 
the M31 disk? To what extent have other recent satellite 
accretions contributed to the growth of the halo? 

As a means to address these questions, we have been 
using the Advanced Camera for Surveys (ACS) on board 
the Hubble Space Telescope (HST) to develop a better 
understanding of the origin and n ature of the stellar sub - 
structure observed around M31. IFerguson et alJ (2005) 
presented the first detailed color magnitude diagrams 
(CMDs) of six regions of prominent inner halo substruc- 
ture, reaching several magnitudes below the red clump 
(RC). Analysis revealed distinct variations in the stel- 
lar populations throughout many of the regions studied, 
consistent with the substructure having more than one 
origin. Two of the regions studied, the NE Shelf and the 
GS, were shown to have remarkably similar populations 
but different line-of-sight distances, con sistent with the 
shelf being a forward wrap of the stream ([Ferguson et al.l 
2005) . A detailed examination of the stellar populations 
in one of the prominent maj or axis overdensit ies, the 
Gl Clump, was presented by iFaria et alJ (|2007l) . This 
field, which lies at ~ 30 kpc, was shown to have expe- 
rienced continuous yet declining star formation during 
the last 10 Gyr an d have a relatively high metallicity of 
[M/H]= -0.4 dex. IFaria et all (|2007l ) argued that these 
properties are more consistent with those of the M31 
outer disk than a low mass accreted dwarf and suggested 
this piece of substructure could have been r ecently torn 
off from the main disk, llbata et al.l (|2005l ) have stud- 
ied the kinematics of many of these features, finding an 
overall pattern of significant rotation. They hypothesize 
that this extended and highly structured rotating struc- 
ture has been formed by accretion. 

In parallel to our efforts to understand the large- 
scale structure and stellar content of M31's outskirts, a 
number of ultra-deep penc il beam pointings have be en 
obtai ned with HST/ACS ([Brown et alJ 120031 [200630, 
2007). With upwards of 32 orbits per pointing, these ob- 
servations have been sufficiently deep to enable stars to 
be resolved to and below the oldest main-sequence turn- 
off, allowing the complete reconstruction of the star for- 
mation history (SFH) in selected locations. iBrown et al.l 
(2006a) present results from inner halo fields they 
associate with the 'Tidal Stream', 'Outer Disk' and 
'Spheroid'. In all cases, the fields were shown to have 
experienced an extended SFH, and differences between 
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39-53-30 


2003-01-08 

iUUJ LI 1 uo 


2430 


5150 




G09458 


A. Ferguson 


00:45:05.0 


39:48:00.0 


2002-10-17 


2430 


5150 


N Spur 


G09458 


A. Ferguson 


00:46:10.0 


43:02:00.0 


2003-07-08 


2475 


5290 


Gl Clump 


G09458 


A. Ferguson 


00:35:28.0 


39:36:19.1 


2003-01-17 


2430 


5150 


Warp 


G09458 


A. Ferguson 


00:38:05.1 


39:37:54.9 


2003-06-10/11 


2620 


5240 


Claw 


GO10128 


A. Ferguson 


00:35:00.3 


40:17:37.3 


2005-01-06 


2469 


5210 


NE Clump 


GO10128 


A. Ferguson 


00:51:56.7 


44:06:38.5 


2004-10-12 


2469 


5210 




GO10128 


A. Ferguson 


00:50:55.2 


44:10:00.4 


2004-10-09 


2469 


5210 


EClJicld 


GO10394 


N. Tanvir 


00:38:19.5 


41:47:15.4 


2005-07-10 


1809 


3000 


GC6.field 


GO10394 


N. Tanvir 


00:38:04.6 


40:44:39.8 


2005-07-23 


1809 


3000 


Brown-stream 


GO10265 


T. Brown 


00:44:18.0 


39:47:36.0 


2004-09-03/05 


2460 


5200 


Brown-spheroid 


G09453 


T. Brown 


00:46:08.1 


40:42:36.4 


2002-12-02/03 


2460 


5200 


Brown-odisk 


GO10265 


T. Brown 


00:49:08.5 


42:44:57.0 


2004-12-20/22 


2460 


5200 



a Total exposure time in F606W filter. Note that the Warp and Brown fields analyzed here do not reach their full depths. 
b Total exposure time in F814W filter. Note that the Warp and Brown fields analyzed here do not reach their full depths. 



them were quantified. Their Spheroid field, which lies ~ 
11 kpc along the southern minor axis, was shown to be 
metal rich (<[Fe/H]> > -0.6 dex) and with a substantial 
intermediate age component. Interestingly, their Tidal 
Stream fi eld, which lies very c lose to the GS field ob- 
served by lFerguson et al.1 (|2005l ). showed almost identical 
populations (but being ~ 1 Gyr younger). Finally, the 
Outer Disk field was best fit by populations of still higher 
metallicities and ages in the range 4-8 Gyr. The interpre- 
tation of these results has since been complicated by the 
suggestion that much of the inner halo, and in particular 
the southern minor axis, is 'contamin ated' by tidal debris 
torn off from the GS progenitor (e.g.. lGilbert et ai1l2007t 
llbata et al.ll2007t ). Indeed, in this paper we highlight the 
global extent of GS progenitor debris across the face of 
M31. 

The halo of M31 is clearly highly complex making 
it difficult to draw firm conclusions about the galaxy's 
assembly history from the study of only a few small 
fields. For example, it has recently been suggested that 
the studies claiming the stellar halo of M31 is markedly 
different from that of the Milky Way may have drawn 
inaccurate conclusions due to comparing kinematically- 
selected Milky Way halo populations with highly pol- 
luted tracts of the M31 minor axis (|Ibata et al.1 120071 : 
iGilbert et al.1 12007). In this study, we have carried out a 
homogeneous analys is of 14 deep HST/A CS pointings, 
including the three iBrown et al.l (|2006al ) fields, which 
probe various regions in the inner halo (R pr oj ^$ 45 kpc) 
and extended disk of M31. Many of these fields were 
specifically targeted as bright and/ or unusual substruc- 
ture identified in the course of the INT/WFC imaging 
survey, while others were taken from studies designed for 
other purposes. In this paper, we present a compara- 
tive analysis of these fields and identify striking trends 
which indicate the inner halo is dominated by a complex 
mix of giant stream material and disrupted disk mate- 
rial. A future paper will present detailed SFH fits to 
these fields. In Section [5] we describe our observations. 
The photometric reduction and completeness testing are 
summarized in Sections [3] and |4] The analysis and re- 



sults are presented in Section [5] Our interpretation of 
the results is discussed in section [6] and summarized in 
section [7] 

2. OBSERVATIONS 

The primary dataset analyzed here comes from our 
programs to obtain deep imagery of prominent substruc- 
ture in the M31 halo (GO 9458 and GO 10128, PI Fer- 
guson). Ten different regions in the M31 outskirts were 
targeted in the course of these studies and nine of these 
are presented in this paper (see Figure 1). Preliminary 
results on the stellar populations at six locations ('Giant 
Stream' (GS), 'North-East Shelf, 'Minor Axis', 'NGC 
205 Loop', 'North ern Spur' and 'Gl Clu mp') were previ- 
ously presented in lFerguson et al.l (|2005[ K Each field was 
observed in the F606W (broad V) and F814W (broad I) 
filters with the Wide-Field Channel (WFC) of the ACS. 
Most of the fields were observed for one orbit in F606W 
and two orbits in F814W, resulting in accurate photom- 
etry to several magnitudes below the horizontal branch. 
The 'Warp' field was exposed for significantly longer how- 
ever the analysis presented here is based on re-drizzled 
images to match the depth of the other fields. The Minor 
Axis, NE Shelf, GS and NE Clump fields all consist of 
two separate pointings to increase star count statistics 
on otherwise low density stellar fields. 

We also utilize two fields ('ECl_field' and 'GC6_field') 
which were observed as part of our imaging survey of 
newly discovered globular c lusters in the outer M31 ha lo 
(GO' 10394, PI Tanvir, see iMackev eTal] (poMl200l '). 
These fields are somewhat shallower than the primary 
fields discussed above (one orbit in each of F606W and 
F814W) but they still provide good sensitivity to well 
below the horizontal branch. With the globular clusters 
masked out (R ~ 30"), the background field populations 
are very useful probes of diffuse structure in M31's ex- 
tended disk. This is the first time results from the Warp, 
NE Clump, Claw, ECl_field and GC6.field have been 
published. 

The remaining three fields in our study come from 
several of the ultra-deep programs of Brown (GO 9453 
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TABLE 2 

Distance and Reddening 



Field E{B-VY N{HI) h R proj c R disk d 

(10 19 cm" 2 ) (kpc) (kpc) 



a Values from the reddening map of lSchlegel etTaTI (1223) ■ 
b M31 column densities; D. Thilker, private communication. 
c Projected radial distance. 

d De-projected radial distance calculated assuming an inclined 

disk with PA= 38.1° and i = 77.5°. 

e Line-of-sight subject to MW confusion. 

and GC 10265). These pointings are labelled 'NGC224- 
Stream', 'NGC224-Halo', and 'NGC224-Disk' in the HST 
archive, though for clarity here we reassign the names 
'Brown-stream', 'Brown-spheroid', and 'Brown-odisk' re- 
spectively. For each field, we retrieved a set of exposures 
from the HST archive and drizzled them to an equiva- 
lent depth of our primary fields. The full-depth pointings 
have already been tho roughly explored and are t he sub- 
ject of several papers (jBrown et al.ll2003l l2006al lbl). We 
note that although the Brown-stream field is sandwiched 
between our own GS pointings, we have chosen to keep 
it separate rather than co-add it for comparative pur- 
poses. The complete observational information for each 
pointing analyzed here is given in Table [1] 

3. DATA REDUCTION AND PHOTOMETRY 

All of the fields were subject to the same method 
of reduction and analysis. Pipeline-calibrated images 
were retrieved from the HST archive and, for each field, 
those images in a given filter were combined with the 
PyRAF TWEAKSHIFT S and MULTIDRIZZLE tasks 
(jKoekemoer et all [20061) . TWEAKSHIFTS was first 
used to calculate any residual offsets between images af- 
ter applying the WCS information from the input head- 
ers. Such shifts were generally found to be < 0.5 pixel in 
the case of images taken within a given visit and < 1 pixel 
for images taken across different visits. MULTIDRIZZLE 
was then invoked to register and stack images using these 
corrected shifts. The final drizzles were conducted using 
the Lanczos3 kernel with pixfrac and scale set to unity. 

Before photometering the images, masks were con- 
structed on the F814W frames to cover diffraction spikes, 
bright background galaxies and saturated stars and then 
applied to the F606W images. Photometry was per- 
formed on the co-added images using the stand-alone 
version of DAOPHOT-II (|Stetsonlll987f ). The stellar den- 
sity in our fields was typically low to moderate. Tests 



revealed that tighter CMD sequences generally resulted 
from aperture photometry alone however full PSF-fitting 
photometry was also carried out as a me ans to discrimi- 
nate b e tween stars and ga laxies (see also iFerguson et aLI 
(|2005f) : iFaria et all (|2007D V After a first pass of object 
detection and aperture photometry, ~250 candidate PSF 
stars were selected from the star catalogue on the basis 
of their isolation, brightness and lack of obvious nearby 
artifacts. In creating an empirical PSF for each field in 
each pass band, any star whose fit was flagged as lying la 
or more from the mean \ 2 value was removed from the 
PSF candidate star list and the PSF re-defined. The ini- 
tial guess (non-varying) PSF model was fit to all sources 
in the catalogue with ALLSTAR II. 

Subsequently, all stars, except PSF candidate stars, 
were fit and subtracted from the image and an improved 
PSF model was built and then re-fit to the entire star 
catalogue. This process was iterated, allowing the PSF 
model to vary linearly with position across the frame, 
until the x 2 statistic of PSF candidate star fits was mini- 
mized and the number of stars retained by ALLSTAR 
II had converged. A further pass of object detection 
and PSF fitting photometry was carried out on a star- 
subtracted frame to reveal stars previously missed. Fi- 
nally, the best-fit PSF model was fit to all stars in the 
catalogue. The photometry list was pruned by rejecting 
stars with any of x 2 , magnitude error or sharpness ly- 
ing outside 3cr of the average value at that magnitude, 
yielding a final list of stellar sources for which we sub- 
sequently adopted their aperture photometry. The two 
catalogues (F814W and F606W) for each pointing were 
then matched to within two pixels. In total, approxi- 
mately 675,000 stars were retained for analysis. 

We then corrected the photometry for aperture effects. 
The 2 to 4 pixel (0."05 to 0."10) radii aperture correc- 
tions were calculated using the average value from the 
photometry of the pruned, matched list of PSF candi- 
date stars. The 4 to 'infinite' pixel aperture correction 
was taken from ISirianni et all ^2005) . Magnitudes were 
placed on the VEGAma g syste m utilizing the zero-point 
values of ISirianni et alJ (|2005[ ) and corrected for fore- 
ground reddening b y interpolating within the maps of 
ISchlegel et all (|1998| ) (see Table EJ. The photometry has 
not been corrected for reddening internal to M31 how- 
ever this correction is not expected to be significant for 
most of our fields since they lie outside the main gas disk. 
Table[2]lists the reddening, internal HI column densities, 
projected and de-projected radial distances of each field 
and is ordered according to increasing projected radial 
distance from the center of M31. We did not correct for 
charge transfer efficiency (CTE) because these data were 
obtained early in the lifetime of the ACS when the cumu- 
lative damage cause d by radiation was minimal (see also 
iBrown et all (|2006af l). Hereafter, all magnitudes refer to 
de-reddened magnitudes in the VEGAmag system. 

4. COMPLETENESS 

The morphology of a CMD reflects the fraction of stars 
in different evolutionary phases; the precise location of a 
star within a given phase is a function of both age and 
metallicity. It is necessary to have a good understanding 
of how complete the CMD is as a function of magnitude 
and color before detailed interpretation. An artificial star 
test algorithm was employed to empirically calculate the 
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Fig. 2. — Results of the artificial star tests on representative fields 
of high, low and intermediate stellar crowding as demonstrated by 
the variation in completeness as a function of magnitude for the 
Spur, Minor Axis and Gl Clump respectively. Also shown is the 
completeness function for the EC1 field which, along with GC6, 
has a reduced exposure time relative to the other fields analyzed 
here. 

photometric scatter and completeness of the data. Each 
image had several thousand artificial stars of a specific 
magnitude (spanning the range 30.0 to 20.0 mag in 0.5 
mag steps) added randomly, where each artificial star is 
a scaled version of the best-fit PSF model. The aim is 
to add as many artificial stars as possible to gain suffi- 
cient statistics, but not so many as to significantly crowd 
the field. It was found that adding ~ 10% of the total 
number of real stars found per band satisfied this condi- 
tion in all fields. The artificial stars are recovered, with 
no prior knowledge regarding their positions, using the 
identical method employed for the real data. They are 
pruned using the same parameters derived from pruning 
the original photometry, not least since the excess ar- 
tificial stars at the magnitude being tested could skew 
the 3<t levels used in the pruning process. The resulting 
catalogue is compared with the actual input positions of 
the artificial stars in order to derive the success rate of 
recovery. Steps are taken to avoid the accidental inclu- 
sion of real stars close to, or superimposed on, artificial 
star positions. Artificial stars coincident with a brighter 
real star are rejected, while artificial stars coincident with 
fainter real stars are retained but considered to have con- 
taminated photometry. The process is repeated until the 
entire magnitude range has been covered. 

Figure [2] shows the completeness rate as a function of 
magnitude for fields spanning a range of stellar densi- 
ties, as well as for one of the shallower fields. For the 
bulk of our fields, completeness levels exceed 80% at 
mF606W = 28.0 mag and tufsuw = 27.5 mag. For the 
slightly shallower EClJield and GC6_field, these limits 
drop to 27.25 mag and 27.0 mag respectively. Thus, a 
high completeness is maintained to generally > 2 mag- 
nitudes below the RC. As expected, the higher density 
fields have somewhat lower completeness at a given mag- 
nitude. Analysis of the photometric errors as a function 
of magnitude and color as derived from artificial star 
tests demonstrate that magnitudes are accurate to <0.10 
magnitudes at tufsuw — 26.5 mag. 



5. ANALYSIS AND RESULTS 
5.1. Color-Magnitude Diagrams 

Figure [3] shows the CMDs of all of the fields rep- 
resented as Hess diagrams with a square-root stretch. 
The total number of stars detected in each field is 
given below the name, and error bars are plotted as 
derived from artificial star tests. The ridge-line of 47 
Tuc (NGC104: [Fe /H] = -0.7, 12.5 Gyr) as taken from 
iBrown et all {2005) is superimposed on the CMDs. It has 
been transferred into the VEGAmag ph otometric sys- 
tem u sing the appropriate zero-points from lSirianni et al.l 
(2005), corr ected to the M31 distanc e modulus of (m- 
M) = 24.47 (jMcConnachie et al.ll200l and de-reddened. 
The CMDs appear broadly similar; all have wide RGBs 
and prominent R Cs typical of M31 popu l ations seen at 
smaller radii (e.g., [M ould & Kristian 1986; Hollan d'et alJ 
1996b iFereuson k, Johnson! 120011: iBellazzini et all 120031 : 
Ferguson et alJ l2005t IBrown et all l2006aD . The RGB 
widths are considerably larger than the photometric er- 
ror bars for magnitudes brighter than mpsuw < 26 mag 
suggesting an intrinsic spread in metallicity. The high 
proportion of RGB stars with colors redder than the 47 
Tuc ridge-line suggests that these stellar populations ex- 
tend to high metallicities ([Fe/H] > -0.7 dex). The Warp, 
GC6_field and Brown-odisk fields reach further to the red 
than other fields. Since these fields are amongst those 
with the highest HI column densities (see Tabled]), it is 
possible that internal extinction is also contributing to 
the color of the RGB. 

On closer inspection, several morphological differences 
are apparent between the CMDs in Figure [31 The most 
striking of these is the morphology of the RC and hor- 
izontal branch (HB) features produced by core Helium 
burning stars. In many of the CMDs, the RC slants 
blueward at fainter magnitudes (hereafter referred to as 
the blue RC 9 ) and is accompanied by an extended hor- 
izontal branch. In contrast, many of the other CMDs 
show fairly round RCs which have a significant luminos- 
ity but little color spread and no evidence for an extended 
horizontal branch. Additionally, these fields also show 
evidence for blue plumes (BPs) of upper main-sequence 
(i.e young) stars. If extended blue HBs are also present 
in these fields, they would partly overlie the BPs mak- 
ing direct detection difficult. One further characteristic 
common to many of the latter set of fields is a promi- 
nent overdensity to the left of the RGB at mpsiiW ~ 
26.1 mag. This feature has the correct magnitude and 
color to be be a main sequen ce turn-off of a 2-3 Gyr 
population (jGirardi et al.ll2000f ). For reasons which will 
become apparent later on, we assign the terms 'stream- 
like' and 'disk-like' respectively to these different CMD 
behaviors. Note that while CMDs within a given group 
share a strong morphological resemblance, they need not 
be composed of exactly identical populations. Indeed, 
the average best-fit age and metallicity of the Brown- 
stream and Brown-spheroid fields, both classified here as 
'stream-like', have b een found to differ b y 0.1 dex and 
0.9 Gyr respectively ([Brown "et~alll2006ah . 

Well populated fields in both categories exhibit a 
prominent RGB bump located below the RC {mpsiAW ~ 
24.6 mag, niF606W — tufsiaw ~ 0.8 mag) implying a 

9 Often referred to as a red horizontal branch in the literature. 
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Fig. 3. — Color magnitude diagrams plotted as Hess diagrams with a square-root stretch to bring out fainter features such as RGB 
and AGB bumps. The bins are 0.025 mag in color and 0.062 mag in magnitu de. The number of stars in each diagram is given below the 
name. The ridge line of 47 Tuc, which has [Fe/H] = -0.7 and age = 12.5 Gyr (Brown ct al. 2005), has been shifted to the distance of M31 
and over-plotted. Photometric errors derived from artificial star tests are shown on the right-hand side of each CMD. The CMDs are > 
90 % complete at mpgi^yy = 27.0mag (~ 80 % for ECl_field and GC6_field). Fields are labelled as disk-like (DL), stream-like (SL) or 
composite (C) as described in the text (Sec 5.2). 



metal rich population. These features are the result of 
stars zig-zagging past the same luminosity value when 
the H burning shell passes the chemical discontinuity left 
behind by the expanding convective envelope causing a 
temporary drop in nuclear efficiency. The RGB bump 
luminosity is expected to decrease with respect to the 
RC as the metallicity increases though it can be difficult 
to interpret in the case of composite stellar populations. 
The most populated fields (NE Shelf, N Spur GC6_field 
and Warp) also feature asymptotic giant branch (AGB) 
bumps at tufsuw ~ 23.1 mag, rriF606W — m F 814W ~ 0.9 



mag. For fields with both RGB and AGB bumps, the fact 
that their luminosities straddle the RC luminosity sup- 
plies an additional clue regarding thei r ages and metallic- 
ities. According to the RC models of lAlves fc~S araicdini 
(1999), such behavior indicates that these populations 
are relatively metal rich. Further, for a single stellar 
population of [Fe/H] = -0.7 dex, this behavior implies an 
age of at least ~5 Gyr, however this restriction is relaxed 
for more metal rich populations 

A few fields exhibit CMDs which share characteristics 
of both stream-like and disk-like behavior. For example, 
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Fig. 4. — Comparisons of individual fields to the GS field. In each panel, the normalized Hess diagram of the field has been subtracted by 
the normalized GS Hess diagram to highlight the relative differences in the underlying populations. Red signifies an under-subtraction and 
blue an over-subtraction. The units are arbitrary and meaningful differences are strong, coherent overdensities in red or blue. Subtractions 
involving the Minor Axis, Brown-stream, ECl_field and NE Clump fields are noisy due to small number statistics. 



the Brown-odisk field has a round RC typical of disk-like 
fields, yet it also shows evidence for a low-level blue RC 
and even an extended HB. Another field with some over- 
lap betw een the two groups is the NGC205 Loop. As 
noted in iFerguson et all (|2005h .this field has a blue RC 
but it lacks the strong RGB and AGB bumps seen in 
other stream-like fields. This is not just a statistical ef- 
fect related to the stellar density on the CMD because the 
RGB bump is noticeably stronger in the Brown-stream 
field which has ^11,000 fewer stars. In addition, there is 
a hint of the bright main-sequence turn-off 'hump' ubiq- 



uitous to our disk- like fields. In Figure 3, the Brown- 
spheroid appears very similar to the GS despite being 
presen ted as pure halo in the literature (|Brown et al.l 
l2006a( bl). However, it will be shown in the following sec- 
tions that there are important differences between these 
fields and that its true nature is likely a mixture of both 
GS and halo populations. These fields thus form a third 
category of 'composite fields'. 

5.2. Differencing the CMDs 

Figure 0] shows the result of subtracting the normal- 
ized GS Hess diagram from the other fields, where the 
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Fig. 5. — Gaussian fits (dashed line) to the red clump mpguw luminosity function (left panels) and color distribution (right panels). 
m RC refers to the peak mfgun/ magnitude and crc to the peak mj?go6w — mF8l4W color. Errors in the fits are ± 0.01 in all cases except 
for the Minor Axis which has an error of ±0.02 in both mjic and cue- 



Hess diagrams have been normalized by the total num- 
ber of stars with F814W band completeness over 90% 
(rriFsiiW < 27.0 mag; corresponding to ^80% in the 
shallower GC6_field and ECl_field). As expected, the 
stream-like fields show excellent agreement with the GS 
leaving negligible residual structure, exce pt in the case of 
the N E Shelf where the RC is brighter. [Ferguson et al.l 
(2005) have argued that this offset in RC magnitude is 
the effect of different line-of-sight distances; we will re- 
turn to this point in Section 16.11 Using their u ltra-deep 
CMD of the giant stream. [Brown et al.l (|2006af ) have de- 
rived an average age of 8.8 Gyr and average metallicity 
of [Fe/H] = -0.7 dex in their field. It is therefore likely 
that this type of population dominates all the stream-like 
fields. 

The subtractions reveal very clear differences between 
the disk-like fields compared to the GS and simple 
line-of-sight distance shifts cannot bring the two into 
agreement. Disk-like fields display a significant over- 
subtraction of the blue RC (at tufsuw ~ 24.3 mag, 0.4 
< mF6Q6W — < 0.9) and an equally prominent 

under-subtraction of the vertically elongated RC. In ad- 
dition, they reveal varying degrees of under-subtraction 
of the blue main sequence populations (rriF606W — 
m<F&i4W < 0.6 mag, mpsiAw > 25.5 mag) indicating 
that these fields have much greater proportions of stars 



with younger ages than the GS. This residual young pop- 
ulation is especially strong in the Warp field. Several 
disk-like fields (the Claw, NE Clump, Gl Clump and N 
Spur) also show a noticeable over-subtraction redwards 
of their RGBs seen as a thin blue strip running paral- 
lel to the RGB for mpsuw < 25.5 mag suggesting that 
GS RGB stars extend to higher metallicities than these 
disk-like fields (note that this over-subtracted band is too 
wide to be caused by age differences alone). 

The composite fields show differing patterns of residu- 
als, none of which exactly mimic those seen in the disk- 
like fields. Despite the NGC205 Loop CMD showing a 
blue RC, there is a strong over-subtraction of that fea- 
ture in Figure U indicating that old metal-poor stars 
constitute a much smaller percentage of this popula- 
tion than in the GS. Indeed the subtraction has revealed 
that the NGC205 Loop field has some disk-like prop- 
erties; a RC elongated in luminosity, a 'hump' consis- 
tent with bright main-sequence turn-off stars (mpsuw ^ 
25.5 mag, niF&o&w — tufsuw ^ 0.6 mag) and an over- 
subtraction redwards of the RGB. In terms of the sub- 
tracted Hess diagram, the NGC205 Loop shows some- 
what similar residuals to the disk-like Claw field, al- 
though the under-subtracted RC is slightly bluer. The 
GS subtracts rather cleanly from the GC6_field although 
residuals confirm that it has an excess of younger and 
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more metal-rich constituents. The Brown-odisk field ex- 
hibits smaller residuals in the region of the RC/HB than 
disk- like fields, implying the possibility of a stream- like 
contribution to the population in this region. 

The last composite field, the Brown-spheroid, displays 
the most unique subtraction signature. Although parts 
of the CMD are well matched by the GS, the red side of 
the RGB is strongly over-subtracted suggesting that this 
field does not contain such high metallicities. Further- 
more, the pattern of RC residuals reveals an excess HB 
sequence not seen in the GS, or any other fields. In fig- 
ures throughout the paper, we will refer to fields as being 
either DL (disk- like), SL (stream-like) or C (composite) 
based on the Hess diagram subtractions. 

5.3. The Intermediate / Old Populations 

The Hess diagram subtractions have confirmed that 
the dominant difference between stream-like and disk- 
like CMDs are those features which reflect the composi- 
tion of the old and intermediate-age core Helium burning 
components. We proceed here to put these differences on 
a more quantitative ground. Figure [5] shows the rriFsuw 
luminosity function (LF) and color distribution in the re- 
gion of the RC for each field. The region selected for the 
color distribution plots was 23.5 < ?tif814W < 24.75 and 
—0.5 < mpQQQw — m F saw < 1-2 to avoid the inclusion of 
BP stars. The dashed line represents the Gaussian fits to 
th ese distributions, usi n g the fitting function formalism 
of iPaczvnski fc Stanekl (jl998l ). At iTiFsuw ^24.0 mag, 
the typical peak RC magnitude, the data are well above 
90% complete and the photometric errors in magnitude 
(< 0.02 mag) and color (<0.03 mag) are small. 

The RC absolute magnitude and color are a function 
of the age and metallicity of a population. For a given 
metallicity, old stars form a fainter RC than young stars. 
The presence of a bright, elongated RC signifies a pre- 
domi nantly metal-rich interm ediate/young aged popula- 
tion (jGirardi &: Sala ris 2001). For a population of known 
age and metallicity, the RC magnitude also acts as a 
standard candle, a fact we will exploit in Section 16.11 
Figure [5] further indicates that the disk-like fields gen- 
erally have brighter and slightly more elongated RCs 
than the stream- like fields, suggesting that they contain 
younger stellar populations in addition to the intermedi- 
ate/old RGB components. In the absence of distance 
effects, the brighter peak magnitude could be the re- 
sult of a contribution from the more evolved counter- 
parts of the blue plume. That is, stars with the same 
age and metallicity but slightly larger masses than the 
young main sequence turn-off stars which causes them to 
evolve into RGB stars more rapidly. The color distribu- 
tions of the RCs are very similar, regardless of grouping, 
though there is marginal evidence for stream-like fields 
having slightly broader color widths. The peak magni- 
tude and color of the Gl Clump {mpmw — 24.12 mag, 
m F606iv — = 0. 75 mag) measured here are in 

excellent agreement with iFaria et al.l (|2007f ) . 

As a test of the robustness of our RC measurements, as 
well as our photometry in general, we measured the RC 
color and magnitude on each individual CMD in those 
fields where we obtained double pointings (the GS, NE 
Shelf, NE Clump and Minor Axis). The peak RC mag- 
nitudes and colors for the two separate pointings, as 
well as for our GS field and Brown-stream, are gener- 
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Fig. 6. — Comparing the completeness-corrected stream-like and 
disk-like RGB luminosity functions. The luminosity functions have 
been normalized at mpgnw = 23.4 mag, between the RC and AGB 
bump. 



ally in excellent agreement, differing by no more than 
±0.03 mag. This indicates that the RC can be used 
as a distance indicator to an accuracy of ±11 kpc. An 
unusual result was recovered in the case of the Minor 
Axis. In that field, the peak RC magnitudes differed 
by a sizeable 0.10 mag in the two different pointings. 
If the populations in these two pointings are identical, 
this magnitude variation would signify a puzzling ~ 38.5 
kpc line-of-sight difference despite a mere ~ 2 kpc pro- 
jected distance between them. More likely, we suspect 
we have serendipitously discovered some small-scale sub- 
structure in this region with intrinsically different pop- 
ulations. The very low stellar density in the individual 
fields (combined these fields have < 17,000 stars) pre- 
cludes a detailed assessment of the similarity of the two 
CMDs, however no obvious differences are apparent. It is 
interesting, however, that the inner of the two pointings 
contains twice as many stars as the outer. 

Figure [6] compares the completeness-corrected RGB lu- 
minosity functions of stream-like and disk-like fields over 
the full luminosity range. The 'stream-like' curve in- 
cludes the GS, Brown-stream, NE Shelf, Minor Axis and 
ECl_field while the 'disk-like' curve includes the Claw, 
NE Clump, Gl Clump, N Spur and Warp. Compos- 
ite fields are omitted. Individual fields in each group 
have first been shifted so that their peak RC magni- 
tude matches the average of that group before being co- 
added and the LFs have been normalized at tofsi4W 
= 23.4 mag, a point mid-way between the RC and 
AGB bump. This plot further emphasizes the gross 
differences between the two populations, including the 
different RC morphologies, with the disk- like function 
peaking at brighter magnitudes. In addition, stream- 
like populations appear to have a higher proportion of 
RGB bump stars than disk-like populations although the 
AGB bumps are similarly concentrated. For magnitudes 
fainter than the RC, the disk-like LF rises more steeply 
than stream-like LF highlighting their young components 
as bright main sequence turn-off stars merge into the 
RGB. 
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Fig. 7. — Close up view of those fields containing blue plume populations. Isochrones from the Girardi ct al. (2000) library of metallicity 
[Fe/H] = -0.4 dex and ages 230 Myr (green), 630 Myr (yellow), 1.0 Gyr (orange), 1.6 Gyr (red) and 2.0 Gyr (magenta) are over-plotted. 
The density of stars falling between adjacent isochrones gives an estimate of the relative strength of star formation in that age bracket. 

be that denser fields would have correspondingly dense 
BPs which is not observed. We also note that disk-like 
fields consistently have higher HI column densities than 
stream-like fields, further supporting the presence of gen- 
uinely young stars in these regions. 

6. DISCUSSION 

6.1. The Nature of the Stream-Like Fields 

The discovery of the giant stream (Ib ata et al.ll200i[) 
was the first spectacular indication that M31 is still ac- 
creting mass. Fields in our study which directly probe 
this substructure are the GS and Brown-stream fields, 
both lying at ~ 20 kpc projected radial distance from 
M31 (see Figure Q}. In addition, we have shown that the 
NE Shelf, EClJield and Minor Axis fields have CMDs 
which are strikingly similar to the GS fields, indicat- 
ing that these fields are also dominated by tidal de- 
bris stripped from the GS progenitor. Although these 
fields generally have smaller projected radii than disk- 
like fields, their positions correspond to larger distances 
i n the plane of the di sk. 

iFardal et all ()2007l ) (hereafter F07) have recently pre- 
sented an N-body simulation of the accretion of a dwarf 
satellite of mass ~ 10 9 M within a realistic M31 poten- 
tial. Their simulation has been tailored to reproduce the 
giant stream and North-east shelf as observed in Figure 
1. As the satellite falls towards the galaxy, long tidal 
streams are produced both leading and trailing the pro- 
genitor core. They identify the trailing stream as the gi- 
ant stream and suggest that the leading stream may have 
wrapped around the inner galaxy at least twice. Figure[8] 



5.4. The Young Populations 

Figure [3] indicates that the disk-like fields, and several 
of the composite fields have sub-populations of young 
stars. Full CMD-fitting is required to place rigorous 
constraints on the recent star formation history of these 
fields but we gain some quantitative insight by compar- 
ing their BPs with the oretical isochrones in Figure 
The lGirardi et alJ (|2000D isochrones of [Fe/H] = -0.4 dex 
with ages of 230 Myr, 630 Myr, 1.0 Gyr, 1.6 Gyr and 2.0 
Gyr best match the BP magnitude and color distribu- 
tion. Lower metallicity isochrones are too blue to match 
the BP while more metal rich isochrones are generally 
too red except in the cases of the GC6_field and Warp 
where there is some overlap. Several fields show BPs 
which reach mFmaw ~ 23 - 24 mag, indicative of stars 
as young as a couple of hundred megayears old. It is 
particularly remarkable that such a young population is 
seen in the NE Clump, a field which sits at ~ 44 kpc (or 
Rdisk ~ 56 kpc ). 

In field populations, blue stragglers are thought to 
form by mass tra nsfer between primordial binaries 
(jDavies et al J 12004). These old blue straggler stars can 
masquerade as you ng main sequence stars on a CMD 
(|Carnev et alJ2005h so it is important to address whether 
the BPs we see c ould be caused by such stars instead. 
iFaria et all (|2007f l provide a thorough discussion of this 
issue for the specific case of the Gl Clump and find that 
while this scenario cannot be completely ruled out, it 
is highly unlikely; these same arguments apply to all the 
fields presented here. If the BPs in our CMDs were made 
entirely of blue stragglers, the naive expectation would 



The Nature Of M31 Halo Substructure 



11 



^ r-t aKsffl; i^n. ' :■ ■< :V:, 




, ± DlHH-lil;* Held 

Stream— lite field 
. Campcnite field 
i I 



TABLE 3 

Line-of-Sight Distances to Stream-like Fields 



£ (dsgre&a} 

Fig. 8. — Locations of our fields with respect to the Fardal et al. 
(2007) simulation showing the predicted distribution of GS progen- 
itor tidal debris around M31. Color coding is as follows: green par- 
ticles represent the GS falling into its first pericentric passage; red 
particles show the fanning out of material into the Northeast shelf 
after the first close pass and approaching their second pericentric 
passage; magenta particles show the fanning out of material into 
the Western shelf after the second close pass and approaching their 
third pericentric passage and blue particles are heading for their 
forth pericentric pass. The locations of our 14 HST/ACS fields 
are over-plotted (see Figure 1 for labels). The ellipse shows the 
approximate extent of M31's bright disk (R = 27 kpc, PA= 38.1° 
and i = 77.5°). The overlap between stream-like fields (open cir- 
cles) and predicted coverage of stream debris is striking. Disk-like 
fields are represented by rilled triangles and composite fields by 
filled stars. 

shows the predicted distribution of the GS progenitor de- 
bris in their model with our HST/ACS pointings over- 
laid. Note that only satellite particles are represented; in- 
side the black ellipse (which indicates the extent of M31's 
bright disk) a significant contribution of M31 disk stars 
is expected. As can be seen, there is an excellent over- 
all agreement between the stream-like fields identified in 
this paper and the predicted pattern of debris. In partic- 
ular, this includes the "Western Shelf" overdensity which 
is probed for the first time by our ECl_field. Figure [8] 
also indicates that much of the minor axis of M31 should 
be contaminated by diffuse GS d ebris stripped off d uring 
various pericentric passages fsee lGilbert et al.ll2007[ ) and 
agrees with our finding that both the Minor Axis and 
Brown-spheroid fields have significant stream-like com- 
ponents. It is also encouraging that our disk-like fields, 
with the exception of the Claw field which projects on 
the very edge of Fardal's Western Shelf, arc found to lie 
well away from regions that contain stream debris in the 
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simulation. 

The three-dimensional distribution of stream debris 
is of great importance fo r placing constraint s on t he 
progenitor orbit (e. g. IMcConnachie et ah ||2003fl ). 
IFerguson et~afl (|2005t ) and lBrown et alj (|2006bh have ex- 
ploited the similarity of the populations in stream-like 
fields to determine different ial distances ba s ed on RC 
magnitudes. In particular, IFerguson et al.l (|2005l ) ap- 
plied this to the GS and NE Shelf fields and derived 
that the latter lay closer to us than the former by a 
factor of 1.07 ± 0.01. The observed RC magnitude is 
distance dependent, so if we assume that the underly- 
ing populations in stream-like fields are the same, we 
can examine their differential line-of-sight distances by 
comparing their peak RC magnitudes. Table [3] presents 
a summary of the distance measurements to stream- 
like fields identified here. The observed distances have 
been determined by assuming that the G S field lies at 
D los = 830 ± 20 kpc (IFerguson et al.ll2005l) with respe ct 
to M31 (Di M = 785±25 kpc: IMcConnachie et all (|2005ft ). 
The RC magnitude differences between the GS and other 
fields then yield relative distances from the center of M31. 
The line-of-sight depth of simulated particles is taken 
from Figure 2b of F07, set on an absolute scale using the 
same Di os (M31). There is an excellent agreement in both 
the differential and absolute line-of-sight distances. 

In using RC brightnesses to calculate approximate line- 
of-sight differences, we have implicitly assumed that all 
the stream-like fields have very similar underlying popu- 
lations. To test the accuracy of this assumption, we have 
used the mean magnitude and color of the RC to shift all 
of the stream-like fields to the same line-of-sight distance 
and reddening as the GS field and then re-subtracted the 
GS Hess diagram (see Figure [9]). These subtractions are 
much cleaner than those in Figure[4]and the residuals are 
consistent with Poisson noise. This lends support to the 
idea that they are composed of the same mix of stellar 
populations as the GS field. The Minor Axis retains a 
slight over-subtraction in the RC area, but it is unclear 
how significant this is given the sparsity of this field. 
Moreover, as we saw in § 15. 3| there may be small-scale 
population variations within this field. 

6.2. The Nature of the Disk-Like Fields 

Aside from sharing a common RC morphology, the 
disk-like fields are also distinguished by the presence 
of a young (< 2 Gyr) population, the prominence of 
which varies from field to field. The NE Clump, Gl 
Clump, Spur, GC6_field and Warp appear to have formed 
stars as recently as a few hundred Myr ago. There is 
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Fig. 9. — Same as Figure[4] but with stream-like fields shifted to 
the same line-of-sight distance as the GS field using RC color and 
magnitude information. 

no clear correlation apparent between the strength of 
the BP population and either the projected or depro- 
jected radial distance of the field, perhaps unsurprising 
given that the gas dis k in M31 is known to be non- 
planar at large radius ([Brinks fc Burton! 119841 ) . How- 
ever, there is a clear correlation between the presence 
of BP stars and the local HI column. For example, the 
fields with the highest unconfused HI column densities 
are the Warp (N H i = 9.8 x 10 19 cm" 2 ) and the Gl 
Clump (Nhi = 26.0 x 10 19 cm~ 2 ) which also exhibit 
the strongest BPs. This could suggest that the young 
stars seen have formed in situ. 

It should be noted that disk-like fields are roughly 
aligned with the major axis and many (Gl Clump, 
N Spur, Warp, Claw and the NE Clump, in addition 
to the NGC205 Loop) probe well-defined high surface 
brightness substru cture lying within the extended disk of 
llbata et all (|2005l ). This intri guing structure is charac- 
terized by strong rotation in the distance range 15 kpc < 
R < 40 kpc, slightly lagging the thin disk component 
by ~ 40 km s" 1 , and has a mode st velocity dispersion 
of 30 km s _1 . llbata et al.l (|2005f ) have discussed vari- 
ous scenarios for the formation of the extended disk and 
concluded the most likely one is via accretion. We revisit 
this important question here in light of the new stellar 
populations constraints we have derived. 

6.2.1. Could disk-like material originate from one or more 
accreted satellite galaxies? 

The idea that the disk-like material results from the re- 
cent accretion of many small satellites provides a ready 



explanation for the jumbled morphology of the substruc- 
ture. However, the homogeneity of the populations in the 
various fields - which are widely separated in distance in 
some cases - is problematic. Within the Local Group low 
mass satellite galaxies generally obey a mass-metallicity 
relation and exhibit a rang e of star formation a nd chemi- 
cal evolution histories (e.g.. lDolphin et aT1l2005l ). Models 
which explicitly take this into account predict that stel- 
lar halos should possess significant age and metallicity 
inhomogeneities resulting from accretion events w ith a 
range of satellite masses ([Font et al.l l2006aL I2007D . As 
we have discussed, such differences a re not obse r ved in 
the disk-like substructure. Moreover, I Font et all (|2007l ) 
have shown that young populations are not expected in 
tidal debris lying within the inner (< 50 kpc) regions of 
hierarchically built halos. 

An alternative scenario is one in which the disk-like 
material comes from a si ngle accretion event inv olving 
a more massive satellite. iPenarrubia et all (|2006h have 
shown that an extended exponential disk structure can 
result from the accretion of a 10 9 — 10 10 M Q object on 
a coplanar prograde object. However, in order to repro- 
duce the high degree of rotation observed in M31's ex- 
tended disk, they require the progenitor to start on a cir- 
cular orbit which is rather at ypical for satellites observe d 
in cosmological simulations (jvan den Bosch et al.lll999h . 
Although the homogeneity of the underlying stellar pop- 
ulation is not a problem in this scenario, the presence 
of genuinely young stars is rather puzzling. Indeed, this 
would seem to require that the disrupting satellite retain 
a substantial reservoir of gas until well within the poten- 
tial of M31, even while much of its stellar component has 
already been stripped. 

Could the source of disk-like substructure be the GS 
progenitor itself? As shown in Table disk-like fields 
have HI column densities at least an order of magnitude 
larger than stream- like fields. It could be that the young 
populations in these fields have formed in situ and are 
unrelated to the dominant stellar component which has 
come from the stream. Our results discount this pos- 
sibility. We have shown that the old and intermediate 
aged stars in stream-like fields have age and metallicity 
distributions which cannot be reconciled with those in 
disk- like fields, as reflected by their differing RC and HB 
morphologies. Furthermore, stream-like fields appear to 
have been more chemically enriched than disk-like fields 
at least a few gigayears before the present time. Cur- 
rent models predict 'first contact' between the GS pro- 
genit or and M31 aroun d ~ 1 Gyr ago (e.g ., llbata et al.1 
12004 iFardal et alJl2006t iFont et al.ll2006bi ) so any mate- 
rial accreted onto the extended rotating disk would al- 
ready have been more enriched than disk-like fields. 

6.2.2. Could disk-like material originate from the thin disk? 

The previous scenarios addressing the origin of disk- 
like material center upon idea that this material has 
been brought into M31 from another system. However, 
the presence of young stars in many extended disk fields 
begs the question as to whether this material could in- 
stead have formed in the gas-rich thin disk of M31 and 
subsequently been torn off and kicked out. 

The Gl Clump is the only disk- like field to have strin- 
;ent c onstraints placed on its SFH s o far. | Faria et ah 
2001 used STARFISH modelling ([Harris fc Zaritskv 
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UPDID of the CMD to infer a high mean metallicity 
([Fe/H] = -0.4 dex) and a large age spread (~ 10 Gyr). 
Although the bulk population is of intcrmediate-to-old 
a S e (~ 6 Gyr), roughly 10% of its stell ar mass has forme d 
within the last 2 Gyr. As discussed bv lFaria et al.l (|2007h . 
these properties are entirely consistent with those of the 
outer disk of M31 where trace a mounts of recent star 
formation are c ommonly observed (|Bellazzini et alj|2003t 
Williams 2001). In conju nction with its kinematic signa- 
ture of disk-like rotation (|Reitzel et a l. 2004; Ibat a et all 
120051) . iFaria et alJ (|2007l ) have suggested a connection 
between the t hin d isk and the Gl Clump. In addition, 
IZucker et al.l (|2004f ) have suggested a link between the 
thin disk and the NE Clump as one possible explanation 
of their observations. In view of the results presented 
here, we explore whether this interpretation can be ex- 
tended to all of our disk-like fields. 

If the disk-like material did originate in the thin 
disk, we require some mechanism capable of restruc- 
turing it and moving material out to large radii (~ 
44 kpc in the case of the NE Clump) and scale 
heights. Recent work has elucidated how interac- 
tions with sub-halos might heat, thicken and perturb 
stellar disks, providing a method capabl e of produc- 
ing s u ch low-latitude sub st ructure (e.g.. iQuinn et all 
19931: IWalker et al.l Il996t IWazauez fe White! 119991 : 
Gauthier et all 120061: IK azantzid is et al.ll2007h . Most re- 
centlv. IKazantzidis et alJ (|2007l ) have performed hybrid 
cosmological plus numerical simulations of dark matter 
sub-halos bombarding a MW type disk galaxy over an 
8 Gyr period and trace the impact on the host galaxy. 
They cull sub-halo properties such as mass functions and 
orbital parameters directly from cosmological simula- 
tions so that sub-halos are representative of their epoch. 
The cumulative effect of six such bombardments, one 
60% as massive as the host's disk, is to produce sev- 
eral long-lived phenomena in the disk such as a signif- 
icant flare, central bar and coherent substructures akin 
to tidal streams. Most damage to the thin disk is caused 
by the most massive satellite. Significantly, although the 
thin disk is thickened and perturbed by these tidal in- 
teractions, forming many messy morphological features 
comparable to those obs erved in M31 and t he M W, it 
survives intact to z = 0. Kazantzid is et all (|2007l ) only 
track the stellar particles in the original thin disk and not 
those of the interacting satellites. In reality, the content 
of the resultant disk will be a complex mix of stars formed 
in situ as well as stars brought in by the satellite. While 
more work is required to establish the viability of this 
explanation, disk heating seems to explain many of the 
observed properties of disk-like fields (and the extended 
disk in general) at least on a qualitative level. 

Many questions arise from this line of thought: Will 
these discrete substructures ultimately disperse to form 
a thick disk in M31? Could the bulk of the damage to 
M31's disk have been caused by the impact of the GS 
progenitor? It is interesting to note that current simu- 
lations suggest the interaction between the GS progen- 
itor and M31 began as recently as ~ 1 Gyr ago (F07). 
This timescale agrees well with the dissolution timescale 
of discret e substructu r es in the extended disk, as esti- 
mated bv llbata et al.l (|2005h . However, there are also 
suggestions that the GS was accreted 6-7 Gyr ago co- 
inciding with the apparent cessation of star formation 



(|Brown et al.ll2006af) . iFont et al.l (|2007l ) present a com- 
prehensive discussion of these two competing scenarios 
based on the results of their numerical simulations. The 
best estimates of the mass of the GS progenitor (10 9 M Q ) 
place it at the low end of the sub -halo mass range fea- 
tured in IKazantzidis et ail (|2007l ). yet it does have a 
highly radial orbit which was found to enhance the dis- 
ruption to the th i n disk . The simulation explained in 
IKazantzidis et ail (|2007h resulted in a final distribution 
of disk stars best described by a thin + thick disk de- 
composition. There is room for much speculation, but 
we are unable to draw conclusions about this within the 
scope of this study. 

6.3. The Nature of the Composite Fields 

Our analysis has revealed that four fields - the Brown- 
odisk, GC6_field, NGC205 Loop and Brown-spheroid - 
do not fit neatly into either stream-like or disk-like cat- 
egories. Several show traits from both. Can these popu- 
lations be described as a simple combination of disk-like 
and stream-like material, or are they genuinely unique 
populations? 

6.3.1. The GC6-field and Brown-odisk fields 

The CMDs of the Brown-odisk and GC6_field exhibit 
many disk-like properties, including prominent round 
RCs and BPs, reflecting their location at the edge of the 
main M31 disk (Rdisk = 25.7 kpc and 26.8 kpc respec- 
tively). However, these same fields also show evidence of 
an extended HB akin to that seen in stream- like fields. 
Furthermore, differencing their CMDs with the GS re- 
sults in lower amplitude over-subtractions than found in 
disk-like fields (see Figure^]) implying a minority stream- 
like component to the stellar populations in these regions. 
If correct, this would be consistent with Figure [5] which 
indicates that both stream debris and thin disk stars are 
expected in considerable proportions at these locations. 

In order to test this further, a composite model CMD 
has been constructed using the N Spur and GS fields 
and compared to the data in Figure [TO] The ratio of 
disk to stream in the model has been scaled by eye to 
produce the cleanest subtraction. Overall, the residuals 
from this subtraction are of significantly lower ampli- 
tude than those in Figure [4] implying that this combina- 
tion of fields is closer in nature to Brown-odisk and the 
GC6_field than the GS alone. In truth, these models are 
crude, but they serve to illustrate that these fields are 
likely disk-like fields contaminated by debris from GS. 
A more accurate breakdown of their populations will be 
achieved via full fitting of the star formation histories in 
future work. 

6.3.2. The NGC205 Loop field 

The NGC205 Loop was so named because, on pro- 
jected surface density maps, it appears as a tidal loop 
emanating from the dE NGC205 galaxy, which sits a 
mere ~ 15 kpc away. Surface brightness profiles of 
NGC205 have long showed significant isophotal t wisting 
implying tidal strippin g of its material by M31 {Hodge 
Il973t IChoi et al J 120021 ). Kinematics of stars thought to 
be associated with the NGC 205 Loop were presented 
by iMcConnachie et"aTI (|2004f ) who identified two kine- 
matic signatures; a majority component attributable to 
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Fig. 10. — Examining the relative contribution of disk and stream 
material to the composite fields. The left-hand panels show the nor- 
malized CMDs of the composite fields while the middle panels show 
model CMDs created from differing amounts of the Spur (or Claw) 
and GS fields. The right-hand panels show the result of subtracting 
these models from the original Hess diagrams. As before, blue cor- 
responds to an over-subtraction and red to an under-subtraction. 

the M31 disk and a small additional kinematically cold 
component (cr.„ ~ 10 kms -1 ) that they associate with 
the NG C205 Loop. Subseq uent examination of the same 
field by llbata et alj (|2005h yielded no evidence for such 
bimodality however and led them to the conclusion that 
the NGC205 Loop was part of the same vast extended 
rotating disk to which the disk-like fields belong (with 
v iag ~ 14 kms -1 , a v ~ 36 kms -1 ). 

We now turn to address what can be learnt about the 
origin of the NGC205 Loop from its stellar populations. 
The NGC205 Loop, despite its stream-like CMD, shows 
significant residuals in the GS subtraction which indi- 
cate a second stellar component is present at this loca- 
tion (Figure [4]) . Is this additional component tidal de- 



bris fro m NGC205 dwarf? iButler k Martmez-Delgadol 
(2005) have recently examined the stellar populations in 
the outskirts of NGC205 and find a population charac- 
terized by a broad RGB, prominent RC and weak HB. 
They estimate a median metallicity of [Fe/H]> —1.06. 
The NGC205 Loop residuals could be consistent with 
such a population. Indeed, the under-subtracted RC in 
this field is slightly bluer than that seen in the disk-like 
fields, suggesting the excess population is more metal- 
poor than the extended disk. On the other hand, a com- 
posite model made from the Claw and the GS (5:1) pro- 
duces a very good match to the CMD of the NGC205 
Loop in Figure 1101 The apparent similarity of the in- 
trinsic NGC205 dwarf stellar population to that of the 
extended disk make it difficult to draw firm conclusions 
about the origin of the NGC205 Loop. Better kinemati- 
cal constraints and more detailed modelling of NGC205's 
orbital trajectory are required. 

6.3.3. The Brown- spheroid field - a glimpse of M31 's 
underlying halo? 

The Brown-spheroid field is unlike any of the other 
fields presented here. Although the GS subtraction is 
fairly clean, it has unveiled an excess extended blue HB 
population not seen in any other subtraction. This sug- 
gests the presence of a higher fraction of metal-poor 
and/or ancient stars than in the GS material. At R pro j — 
11.5 kpc, this is the innermost of our ACS pointings; if a 
weak smooth halo component exists in M31, then it is in 
t his field that we wou ld expect the strongest signature. 

Gilb ert et alj (|2007T ) have recently revisited the nature 
of the Brown-spheroid field in light of their spectroscopic 
study of stars along the minor axis of M31. Based on 
kinematics, they estimate that ^44% of stars in this field 
belong to a cold component with the same metallicity as 
the GS and they have shown that this fraction agrees well 
with the predictions of the F07 simulation. Figure [10] 
shows the result of subtracting a 44% GS component 
from the Brown-spheroid field. The residual population 
is characterized by a RC which has narrowed somewhat 
in luminosity, the lack of an RGB bump and the presence 
of a prominent blue HB. It is tempting to speculate that 
this represents the true underlying halo of M31. 

Our own Minor Axis field, which sits ~ 8 kpc fur- 
ther out than the Brown-spheroid field, was originally 
intended to probe the uncontaminated halo, however we 
have found here that it too s hows the signature o f con- 
tamination by stream debris. iGilbert et al.l ([2007) trace 
the kinematics in a location close to our Minor Axis field 
finding that ~31% of the stars are part of a cold giant 
stream-like component. It may appear surprising that 
this field compares so much better to the GS than the 
Brown-spheroid when the fraction of stream present is 
estimated to be even lower. However, it should be noted 
that the Minor Axis field lies at a very similar projected 
radial distance from the center of M31 as the GS field 
(see Table 2). Therefore, if the underlying M31 halo is 
spherical, the fraction of ancient metal-poor components 
in these two fields should be almost identical and thus 
cancel out cleanly in subtractions between the two. 
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7. SUMMARY 

We have presented the largest and most detailed sur- 
vey to date of stellar populations in the outskirts of 
M31. We have carried out a homogeneous analysis of 14 
deep HST/ACS pointings spanning the range 11.5 kpc 
<Rp ro j S 45 kpc. These pointings sample well-defined 
substructure identified in the course of the INT/WFC 
imaging survey, as well as the more diffuse extended disk. 
We have shown how stellar populations in these fields can 
be generally divided into two clearly distinct categories 
based on their CMD morphologies: stream-like or disk- 
like. Our analysis lead us to the following conclusions: 

• Stream-like fields have blue RCs and extended hor- 
izontal branches yet are more metal enriched than 
disk-like fields. They show no evidence of re- 
cent star formation. Our analysis reveals for the 
first time that the Western Shelf (probed by the 
ECl_field) is consistent with being material torn off 
from the GS progenitor. We compare the spatial 
and line-of-sight distribution of stream-like fields 
with the F07 simulation of the GS progenitor or- 
bit and find an excellent agreement. In this pic- 
ture, the GS and Brown-stream fields sample the 
giant stream itself while the NE Shelf and EC 1 -field 
probe the Northeastern and Western Shelves. The 
Minor Axis and Brown-spheroid fields are sig- 
nificantly contaminated by stream material that 
has undergone several pericentric pa ssages, as pre- 
vious ly suspected from kinematics (|Gilbert et al.l 
l20f)7h . 

• Disk-like fields have rounder RCs which have sig- 
nificant luminosity width and show no evidence 
of an extended horizontal branch. All show ev- 
idence for young (< 2 Gyr) populations, and in 
some cases star formation as recent as 250 Myr 
ago. Disk-like material could have either an ex- 
ternal or an internal o rigin. Recent work (e.g., 
iKazantzidis et all I2007D describes how such struc- 
tures resembling tidal debris can be created from 
interactions between the stellar thin disk and mas- 
sive satellites on close pericentric orbits. Given 
the uniform populations in these fields, including 
the ubiquitous presence of young populations, and 
the high rotation observed in th e extended disk t o 
which the disk- like fields belong (|Ibata et al.ll2005l) . 
we favor this scenario over the original assertion 
that this vast extended rotating component has 
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formed via accretion. 

• Several fields are identified as being 'composite' in 
the sense of having both stream-like and disk-like 
properties. Amongst these are the Brown-odisk 
and Brown-spheroid fields which have been the tar- 
gets of ultra-deep HST programs aimed at deter- 
mining the SFH o f the M31 disk and h alo back to 
the earliest times (|Brown et~aT][2006ah . Our find- 
ings indicate that the stellar populations in these 
regions are probably much more complicated than 
originally envisioned, making interpretation of the 
resulting SFHs difficult. Another composite field 
is the NGC205 Loop. We have shown that this 
field likely contains GS contamination in addition 
to an underlying component which is consistent 
with being either material from NGC205 and/or 
from M31's extended disk. 

We have shown that much insight can be gained into 
the nature and origin of substructure in M31 from the 
comparative analysis of CMDs reaching a few magni- 
tudes below the RC. As new surveys uncover an ever- 
increasing level of complexity in this system, a shallow- 
depth multi-pointing approach such as this may prove 
very valuable. In forthcoming work, we will present de- 
tailed SFH fits in the various fields analyzed here as a 
means to confirm and better quantify their similarities 
and differences. It will also be of great interest to extend 
this analysis to the t idal streams recen tly discovered in 
the far outer halo bv llbata et al.l (|2007t ). 
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